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Rapid adiabatic preparation of injective PEPS and Gibbs states 
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We propose a quantum algorithm for many-body state preparation. It is especially suited for 
injective PEPS and thermal states of local commuting Hamiltonians on a lattice. We show that for 
a uniform gap and sufficiently smooth paths, an adiabatic runtime and circuit depth of 0(polylog N) 
can be achieved for 0{N) spins. This is an almost exponential improvement over previous bounds. 
The total number of elementary gates scales as 0(A polylog A). This is also faster than the best 
known upper bound of O(A^) on the mixing times of Monte Carlo Markov chain algorithms for 
sampling classical systems in thermal equilibrium. 


Quantum computers are expected to have a deep im¬ 
pact in the simulation of quantum many-body systems, as 
initially envisioned by Feynman [1]. In fact, quantum al¬ 
gorithms have potential applications in diverse branches 
of science, ranging from condensed matter physics, atom 
physics, high-energy physics, to quantum chemistry [2]. 
Lloyd [3] was the first to devise a quantum algorithm to 
simulate the dynamics generated by few-body interact¬ 
ing Hamiltonians. When combined with the adiabatic 
theorem [4, 5], the resulting algorithms allow one to pre¬ 
pare ground states of local Hamiltonians, and thus to 
investigate certain quantum many-body systems at zero 
temperature. Quantum algorithms have also been intro¬ 
duced to prepare so-called projected entangled pair states 
(PEPS) [6-8], which are believed to approximate ground 
states of local gapped Hamiltonians. Furthermore, quan¬ 
tum algorithms have also been proposed to sample from 
Gibbs distributions [9-14], which describe physical sys¬ 
tems in thermal equilibrium. The computational time of 
most of these algorithms is hard to compare with that 
of their classical counterparts, as it depends on specific 
(e.g., spectral) properties of the Hamiltonians which are 
not known beforehand. However, they do not suffer from 
the sign problem [15], which indicates that they could 
provide significant speedups. 

Quantum computers may also offer advantages in the 
simulation of classical many-body systems. For instance, 
quantum annealing algorithms [16, 17] have been devised 
to prepare the lowest energy spin configuration of a few- 
body interacting classical Hamiltonian, which has obvi¬ 
ous applications in optimization problems. Quantum al¬ 
gorithms have also been proposed to sample from their 
Gibbs distributions at finite temperature [18-23]. Apart 
from applications in classical statistical mechanics, simi¬ 
lar problems appear in other areas of intensive research, 
e.g., machine learning. Speedups as a function of spectral 
gaps have been analysed in Refs. [12, 21, 22]; the scaling 
with large system sizes, which is of particular interest for 
applications in deep machine learning [24] , is however not 
optimal. 

In this Letter we propose and analyse a quantum al¬ 
gorithm to efficiently prepare a particular set of states. 
This set contains two classes relevant for lattice prob¬ 


lems: (i) injective PEPS [25]; (ii) Gibbs states of locally 
commuting Hamiltonians. Glass (ii) contains all classical 
Hamiltonians, and thus the quantum algorithm allows 
us to sample Gibbs distributions of classical problems at 
finite temperature. 

Our algorithm outperforms all other currently known 
algorithms for these two problems in the case that the 
minimum gap A occurring in the adiabatic paths (to be 
defined below) is lower bounded by a constant. We show 
that the computational time for a quantum computer, 
given by the number of elementary gates in a quantum 
circuit, scales only as 

r = 0(A polylog (A/e)), (1) 

where A is the number of local Hamiltonian terms, e 
the allowed error in trace distance and the degree of the 
polynomial depends on the geometry of the lattice. Note 
that an obvious lower bound on the computational time 
is H(A), as each of the spins has to be addressed at least 
once. Thus, Eq. (1) is almost optimal. Furthermore, the 
algorithm is parallelisable, so that the depth of the circuit 
becomes 

D = 0 (polylog(A/e)). (2) 

This parallelisation may also become very natural and 
relevant in analog quantum simulation, as is the case for 
atoms in optical lattices [26]. 

One of the best classical algorithms to sample accord¬ 
ing to the Gibbs distribution of a general classical Hamil¬ 
tonian is the well-known Metropolis algorithm [27]. The 
currently best upper bound to its computational time is 
T = O(AVAstoch) [28], where Agtoch is the gap of the 
generator of the stochastic matrix. We will see that given 
any stochastic matrix, one can always construct a quan¬ 
tum adiabatic algorithm with the same gap A = Agtoch, 
and thus we obtain a potential quantum speedup of al¬ 
most a factor of A. Under parallelisation, the circuit 
depth is almost exponentially shorter. Our algorithm 
to prepare injective PEPS also provides a better scaling 
than the one presented in Ref. [7]. 

The class of states we consider in this Letter can be 
thought of as commuting finite range operators acting on 
a set of maximally entangled pair states (Fig. 1). More 
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FIG. 1. The general class of states, Eq. (3). Finite range 
operators (red) acting on a collection of maximally entangled 
pair states (blue) distributed on a graph. 


precisely, consider a regular lattice in some dimension, 
and let Q = {V,£) be the associated (infinite) graph. 
We endow Q with a distance d, the minimum number 
of edges separating two vertices in V. We associate a 
d-dimensional Hilbert space, Hv, to each of the vertices 
V £ V. Consider the set A of interaction supports, i.e., 
A is a collection of sets of vertices whose relative dis¬ 
tance is at most a constant R, the interaction length, 
and consider for each A G A an interaction Qx which 
is an operator supported on We assume that 

they are strictly positive, 1 > Q\ > gol; and mutually 
commute, [Qx,Qx'] = 0. Consider also a set Y of mutu¬ 
ally excluding pairs of neighbouring vertices. Moreover, 
let Ajv be a finite subset of A with |Ajv| = N, and define 

\4'n) oc JI Qa 0 (3) 

where T 7 v = {^GY|/rn (UagAjv ^ 
of pairs with a vertex in A^r, and K*) i® 

an unnormalized maximally entangled state between the 
pairs of vertices in Tjy. We will give a quantum algorithm 
to prepare the state Eq. (3), and analyse the runtime as 
a function of N and other spectral properties. In the 
following, we drop the subindex N to ease the notation. 
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FIG. 2. (a) Projected entangled pair states, (b) Purification 
of a thermal state. For each system qudit, we introduce an 
ancilla to be placed in a maximally entangled pair with its 
system particle, then apply to the system. 


As mentioned above, Eq. (3) includes two relevant 
classes of states. The first is the class of injective PEPS. 
The graph is composed of nodes, each of them includ¬ 
ing a set of vertices (Pig. 2a). In this case, T contains 
pairs of vertices in nearest neighbor nodes, whereas A 
contains each node. The operators Qx act on different 
nodes, and therefore trivially commute. The resulting 
state is just a PEPS, which is injective since each Qx is 
invertible. In fact, every injective PEPS can be expressed 
in this form up to a local unitary using a QR decompo¬ 
sition. The second class is the class of Gibbs states of 
commuting Hamiltonians [29]. To see this, consider the 
graph which contains sites composed of two vertices, one 
of them is called “system” and the other “ancilla”. The 
set T contains all sites, whereas A contains interacting 
system vertices (Fig. 2b). The relation with Gibbs states 
is evident if we write Qx = where ||/ia|| < 1, 

and take into account that they mutually commute. It is 
easy to see that if we trace the ancillas, we obtain 

poce"^^, (4) 


where H = I^agA^a- 

The state Eq. (3) is the unique ground state of a 
frustration-free local Hamiltonian that can be written as 

G = ^ (5) 

/xGT 

with 

gm=( n qaM ( n qaM . ( 6 ) 

where A^ = {A S A | A n p 7 ^ 0} is the set of supports 
whose interactions act nontrivially on p, and Pf^ is the 
projector onto the subspace orthogonal to \4>^)fj.- Notice 
that since each is supported in a region of radius R 
around /i, G is indeed local. 

The state Eq. (3) can be prepared using an adiabatic 
algorithm. For that, we define a path Qx{s) with unique 
ground state |^(s)), where s £ [ 0 , 1 ], with Qa(0) = 1 and 
(5 a(1) = Qx- We can choose (5 a(s) = (1 — s)l + sQx- 
In the case of the thermal state, we can also choose 
Qxis) = Then, by starting with |())(0)) and 

changing the parameter s : 0 —>■ 1 sufficiently slowly, we 
will end up in the desired state |(/)(1)). The runtime for 
this preparation, as measured by the number of elemen¬ 
tary quantum gates, is unpractical, however, as it scales 
a.sT — polylog(A/£)), where e is the toler¬ 

ated error and A is the minimum spectral gap along the 
path. Indeed, the adiabatic theorem [30] gives an adia¬ 
batic runtime of r = so that Hamiltonian 

simulation [31] gives T = 0{tN‘^ po\ylog{N/ e)). 

To obtain a better scaling, we first use a variant of 
the adiabatic theorem with almost exponentially bet¬ 
ter runtime dependence on the error using a sufficiently 
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smooth reparameterisation of the Hamiltonian path. The 
quadratic scaling of the runtime with the derivative of 
the Hamiltonian, however, leads to an unpractical depen¬ 
dence on N since the Hamiltonian contains 0{N) terms 
that change with time. To avoid this, we change the 
Q’s individually, leading to an adiabatic runtime of r = 
0(fVlog^^“ (A^/eA) A“^). This, however, uses Hamilto¬ 
nians acting on the whole system, despite only the change 
of a single Q, which would result in an additional factor 
of polylog(A/e)) for the computational time mea¬ 

sured by the number of elementary gates. We circumvent 
this problem by using Lieb-Robinson bounds [32] and the 
frustration freeness to show that under the assumption 
of a uniformly lower bounded spectral gap, it is at each 
step sufficient to evolve with a Hamiltonian acting only 
on 0(polylog(A/£:)) sites instead of the full lattice. 

Thus, define a sequence of N Hamiltonian paths by 
enumerating the elements of A as Ai,..., Aat, and define 

Gnis)= Y. ( 7 ) 

fj,er 

d(/i,A„)<xlog^+“(Af/£) 

for n = 1,..., A, where x is a constant control parame¬ 
ter, and 


Gn,tJ.{s) 


with 








An,m ( 5 ) 


( 8 ) 

m<n 

WA„(s) m = n (9) 

[(5a„,(0) = 1 m> n. 


Notice that Gn is supported on a region of radius 
0(log^'''“(A/£)) and dG^/ds is supported on a region 
of bounded size. By reparameterising G„(s) —?► G„(/(s)) 
with / a function in the Gevrey class l-|-a, we can assume 
Gn to be in the same Gevrey class [33]. 

Consider the sequence of Schrbdinger equations 


*dt ~ ) I'*/'") > IV'n-n(O)) = \tpn{Tn)) (10) 


for times = 0(log(A/e)^+“), stating in |V’i(0)) = 
|(()(0)), the trivial ground state of Gi(0). The al¬ 
gorithm proceeds by running Hamiltonian simulation 
[31] on this sequence of adiabatic evolutions. Since at 
all times we only evolve with Hamiltonians acting on 
0(polylog(A/e)) sites, the number of gates only grows 
as Eq. (1). Moreover, the evolution of consecutive G„s 
can be parallelised if their support is disjoint, i.e., if 
Gn, • ■ ■, Gn+i have disjoint supports, the subsequence can 
be replaced by their sum without altering the evolution. 
Since |suppG„| = 0(polylog(A/e)), it is clear that an 


ordering of the A„ can be chosen such that subsequences 
of length fl(A/polylog(A/e)) of the G„s can be paral¬ 
lelised at a time, resulting in a circuit of depth Eq. (2), an 
almost exponential improvement over previous bounds. 

In the following, we show that for a uniformly lower 
bounded gap, the error of the above algorithm is bounded 
by e. First, we use that under sufficient smoothness con¬ 
ditions on a Hamiltonian path G(s), the final error can 
be almost exponentially small in the adiabatic runtime. 
Indeed, as proven in the Supplemental Material, if G is 
in the Gevrey class \ + a and d^G/ds^ = 0 at s = 0,1 
for all A: > 1, then an adiabatic runtime of 


T = O 





( 11 ) 


is sufficient for an error e, where A is the minimum gap 
of G(s) and K = \ supp dG/ds| if G is local. The required 
smoothness conditions can always be achieved with a 
suitable reparametrisation of the path G(s) —>■ G(/(s)). 

This allows us to implement the global change of the 
Hamiltonian, Eq. (5), as a sequence of N local changes. 
Define the sequence of Hamiltonian paths. 


Gnis) = Y GnAs)- ( 12 ) 

mgt 


Notice that Eq. (12) is the same as Eq. (7), but con¬ 
tains all local terms Gn,ii- The weak dependence on 
in Eq. (11) ensures that the accumulated error 
under the sequential evolution with Eq. (12) remains 
small. Indeed, for a final error e, it is sufficient that 
the evolution with each G„ in this sequence only gen¬ 
erates an error of at most e/N. Equation (11) and 
|suppdG„/ds| =0(1) imply that this can be achieved 
in a time Tn = O (\og^'^°‘{N/eAn)AA), where A„ is the 
minimum spectral gap of Gn- A decomposition into a 
circuit then requires T = O (A^ polylog(A/£A) A“^) el¬ 
ementary gates, where A = min„A„. This is already 
an improvement by a factor A over the naive change of 
the entire Hamiltonian, assuming similar behaviour of A 
compared to the spectral gap of the original path G(s). 

Assuming that A = D(l), we can further improve on 
this using Lieb-Robinson bounds to localise the effect of 
the adiabatic change. Indeed, we show in the Supple¬ 
mental Material that local terms in Eq. (12) which are 
supported at a distance D(log^'''“(A/£)) away from the 
support of dG„/ds do not significantly contribute to the 
unitary evolution generated by Eq. (12). This allows 
the replacement of Eq. (12) with Eq. (7) without sig¬ 
nificantly altering the evolution and thus the final state. 
Notice that G„ only acts on 0(polylog(A/£)) sites and 
Tn = 0(polylog(A/£)) for all n. Thus, its unitary evolu¬ 
tion can be simulated with only 0(polylog(A/£)) gates. 
Hence, we hnally obtain a number of gates in the circuit 
model that grows only as Eq. (1) for a constant error and 
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lower bounded spectral gap. Using the described paral¬ 
lelisation, we finally obtain a circuit depth Eq. (2), as 
claimed. 

In the analysis above, we have assumed a gap > A 
along all N paths. This assumption can in fact be relaxed 
to a gap at either s = 0 or s = 1 (see Supplemental 
Material), using the positivity condition on Q\. We thus 
say that the system has a uniformly lower bounded gap 
A if for all finite subsets A C A, the Hamiltonian Eq. (5) 
has a spectral gap > A. Under this assumption [34], the 
circuit depth Eq. (2) can be guaranteed [35]. 

For the preparation of thermal states of classical 
Hamiltonians H, it is natural to compare these re¬ 
sults with the performance of classical Monte Carlo 
Markov chain algorithms for Gibbs sampling such as the 
Metropolis algorithm or Glauber dynamics. Notice that 
due to the nature of their implementation, a fair compar¬ 
ison of performance should compare the mixing time of a 
discrete-time Markov chain to the number of elementary 
quantum gates, whereas the mixing time of a continuous¬ 
time Markov chain should be compared to the circuit 
depth. The best known upper bound on the discrete 
mixing time for Monte Garlo Markov chain algorithms 
for sampling from Gibbs distributions of classical Hamil¬ 
tonians given just the promise of a spectral gap scales 
as 0{N^). Although under certain additional assump¬ 
tions such as translational invariance [36, 37], weak mix¬ 
ing in two dimensions [38], or high temperature [39], the 
existence of a logarithmic Sobolev constant and hence 
the rapid (discrete) mixing time of 0{N log N) can be 
proven, no such proof exists for the general case to the 
best of our knowledge. Our scheme thus outperforms 
classical Monte Garlo algorithms whenever rapid mixing 
cannot be shown even in the presence of a uniform gap. 

Note that any classical Monte Garlo algorithm can be 
realised as an adiabatic algorithm, as has, e.g., been ob¬ 
served in Ref. [20]. Indeed, if S is the generator matrix 
of a continuous-time Monte Carlo algorithm that satisfies 
detailed balance with respect to the Gibbs distribution, 
G = —is Hermitian. This Hamiltonian 
has the same spectrum as —S' and has the unique ground 
state J-l-)®’^. For classical Hamiltonians H, this 

state has the same measurement statistics as pp for ob¬ 
servables that are products ol az- By introducing an an- 
cilla for every particle and applying the map ji) i—)■ \ii), 
the purified version of the thermal state can also be recov¬ 
ered, and its parent Hamiltonian has the same spectrum 
as —S within the symmetric subspace. Thus, any classi¬ 
cal system with a uniform spectral gap for the generator 
matrix can be turned into a rapid adiabatic algorithm. 

For quantum Hamiltonians, notice the restriction to 
commuting local terms. For noncommuting local terms, 
an approximate quasilocal parent Hamiltonian can be 
considered above some constant temperature that allows 
the preparation in polynomial time. We describe this 
procedure in the Supplemental Material. 


For the preparation of injective PEPS, the given al¬ 
gorithm is similar to Ref. [7], which, however, requires 
a runtime of O(A^) in the uniformly gapped case, due 
to the use of phase estimation and the “Marriot-Watrous 
trick”, which are computationally expensive for large sys¬ 
tems. The better runtime of the present algorithm is 
largely due to replacing these subroutines by a local adi¬ 
abatic change. 

Throughout the analysis of this Letter, we focused 
on the case where a uniform constantly lower-bounded 
spectral gap is assumed. This assumption is only used 
to obtain a small number of elementary gates and 
circuit depth, whereas the adiabatic runtime of r = 
0(Alog^^“(A/eA)A“^) is independent of this assump¬ 
tion [40]. In contrast, the runtime of the algorithm to 
prepare PEPS given in Ref. [7] only grows as T ^ A“^ for 
small gaps, and for thermal states, algorithms based on 
quantum walks, phase estimation, and the quantum Zeno 
effect have been proposed with a runtime of T ^ A“^/^ 
[12, 21, 22], albeit with worse scaling in the system size. 
We believe that similar techniques can be applied to our 
scheme of local changes to obtain a good scaling of the 
runtime for both large system sizes and small spectral 
gaps. Moreover, it would be interesting to investigate if 
this scheme of local changes can also be applied to speed 
up classical Monte Garlo algorithms. 

We have also shown that the algorithm can be par¬ 
allelised, thus giving rise to a circuit depth that scales 
only polylogarithmically with N. This is particularly at¬ 
tractive for certain experimental realisations of analog 
quantum simulators, such as with atoms in optical lat¬ 
tices [41] or trapped ions [42], where this could be carried 
out in a very natural way. 

We thank A. Lucia, D. Perez-Garcia, A. Sinclair, and 
D. Stilck Franga for helpful discussions. This work was 
supported by the EU Integrated Project SIQS. 
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SUPPLEMENTAL MATERIAL 


Proof of the adiabatic theorem with almost exponential error decay 


In this section, we prove a variant of the adiabatic theorem that only requires a runtime almost exponentially small 
in the allowed error. Our proof largely follows the proof given in [45], which is based on the method of adiabatic 
expansion [46]. The adiabatic expansion in [46] establishes an approximation of the time-dependent Schrodinger 
evolution in terms of the instantaneous ground state and its derivatives, but on its own does not necessarily imply an 
adiabatic theorem because it assumes a special initial state. Our proof, like [45], resolves this problem by exploiting 
the Gevrey-class condition which allows to satisfy these initial conditions, and uses this expansion to establish a bound 
on the required runtime. However, unlike [45], which only proves the almost exponential dependence of the runtime 
with respect to accuracy, our proof also explicitly establishes the dependence on all other parameters such as the 
spectral gap and the bound on the Hamiltonian derivatives [47]. 

Consider a smooth path of Hamiltonians, G(s), s G [0,1], acting on a finite-dimensional Hilbert space "H. Let (/>(s) 
[48] be the ground state of G(s) and tj}{e, s) the solution of the following Schrodinger equation: 

*eV’(e,s) = G(s)?/'(e,s), V'(O) = (13) 


where 1/e = r is the runtime of the adiabatic algorithm, and ' denotes derivative with respect to s. We assume 
furthermore that the ground state energy of G(s) is 0 (i.e., we fix the phase of ijj) and that it has a gap at least A 
throughout the whole path. By an appropriate choice of the phase of </, we can without loss of generality assume that 
0(s) ) = 0. In the following, we will sometimes drop the explicit dependence on s to simplify the notation. 


Unless otherwise stated, ]].]] will always denote the operator norm for operators and the Euclidean vector norm for 
vectors (it will always be clear from the context which one is used). In this section, we prove the following theorem. 


Theorem 1. Suppose that all derivatives of G vanish at 0 and at 1, and moreover that it satisfies the following Gevrey 
condition: there exist non-negative constants K, c and a such that for all k > 1, 


Then, 


| 1 G (*^)|1 < Kc '^ 


[fc!]i+“ 
(fc + l)2' 


mm 1) - e*V(l)|| < 8ce^ 




(14) 


(15) 


Notice that we don’t require the Gevrey condition (14) to hold for k = 0. Therefore, in the application of Theorem 1 
in the main text, K = 0(1) since along the paths G„(s) (as defined in (12) in the main text), only 0(1) local terms 
change. 

Proof of Theorem 1. Following the adiabatic expansion method from [45, 46], we search 'tf{e,s) in the form of an 
asymptotic series expansion by constructing vectors </j(s), s G [0,1], j > 0, such that for all M > 0, 


M-l 

3=0 




= 0(e"). 


(16) 


We first show an explicit expression for (fj provided that such an expansion exists. Second,we prove that the expansion 
really exists if G^^^ (0) = 0 for all k by giving an explicit error bound. Third, to connect the expansion to the adiabatic 
theorem, we show that 


M-l 

mm ]]?/(£, 1) - e*V(l)II <2||^/(e,l)- ^ <(>i(s)e^||, (17) 

3=0 

for some 6 for all M if G*^*^(l) = 0 for all k. This already proves an error bound of 0{e^) for any M. Finally, if 
G is Gevrey class, then the error bound can be expressed with the help of the parameters appearing in the Gevrey 
condition and using a suitable choice of M yields to the bound in Eq. (15). 
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Explicit form of (pj. To satisfy the equation at s = 0, we require (^o(O) = ^(0) and 0j(O) = 0 for all j > 0. 
Furthermore, substituting back the ansatz to the Schrodinger equation Eq.(13), following [46], we arrive at the 
recursion 


= tpj(j) + iG ipj=i / dt'G , 


for all j > 0, where ^j{s) is a complex number and G ^ is the pseudo-inverse of G, and initial values are 

</>o(s) = (pis) 

(^o(s) = 1- 


(18) 


(19) 

( 20 ) 


Note that (/5j(0) has to be zero in order for (t>j{Q) to be zero, but this is not a sufficient condition. We will investigate 
below when (j)j{0) = 0 can be satisfied. 

Existence of the expansion. To satisfy (^j(O) = 0 for j > 0, (pj-iiO) needs to be parallel to (^(0). This is satisfied if 
all derivatives of G are 0 at s = 0 (see Lemma 2 below). We show that if this condition is fulfilled, then the expansion 
exists. 

Define the truncation of the asymptotic series expansion, 


M-l 

'ipM= ^3^^- 

3=0 


( 21 ) 


Note that if II?/; —then the expansion exists. Indeed, then ||?/; —?/;m|| = ||V^ —^/’M-i-i+e^</>M|| = O(e^). 
By construction, ipM almost satisfies the Schrodinger equation: ieipM = GipM +ie^vM-i- Let U{s) be the dynamics 
generated by G/e. Then, ||?/;M(e, s) — ^/^(e, s)|| = ||t/(s)VM(e, s) — (/)(0)|| and 


\U{s)^'iljM{e,s) - (/i(0)|| = 




M-l 


ds'U^M-i 


< e 


M-l 


ds' 


0M-] 


( 22 ) 


where we used that if the first M derivatives of G are 0, then = 4>{0)- This proves the existence of the 

expansion. 

Connecting the expansion to the adiabatic theorem. Using the triangle inequality, we obtain 

min||?/;(e, 1) - e*®(/)(l)l| < ||V’(e, 1) “'0M(e, 1)|| -f min ||?/;M(e, 1) - e*V(l)ll • (23) 

6 6 


In Lemma 2, we prove that if the first M derivatives of G(s) vanish at s = 1, then (/)j(l) is parallel to </>(!) for 
all J = 1,..., M. Therefore, '0M(e, 1) is parallel to ^i>(l), so that mine ||V'M(e, 1) — e®®^(l)|| = |||?/;M(e, 1)|| — 1|- But 
1 = \\'tp{e, 1)11, so using the triangle inequality, we get mine ||l/’M(e, 1) — e*®(/>(l)|| < 1) — V'(f) 1)11- Therefore, 

min||?/;(e,l)-e*V(l)|| < 2||?/;(e, 1) - ?/;M(e, 1)||. (24) 

u 


Choice of M. From Eq. (22) and (24), 


min ||V;(e, 1) - e*V(l)|| < 2||?/;(e, 1) - V'M(e, 1)|| < 2e 


M-l 


ll0M-l||, 


(25) 


so we only need to bound ||(/)m-iI 1- We do this by using that G is Gevrey class. From Lemma 7 below, 

M-l 


||0M-i|| < 2-^ • 2c— ( 


A3 I 3 


[M!] 


1+a 


{M + lf 


(26) 


Therefore, 


M-l 




A3 ^ 3 y (M+l)2- 


IIV'(e,l) - <(’(1)11 < £■ 


(27) 






















Changing M to M + 1 and using that [{M + l)!]^+“/(M + 2)^ < we obtain 


K /dTT^^ ^ 

||V>(e,l)-<^(l)||<8c-( — 


K"^ , 2 /^47r2V 


M 


This is true for any M, so setting 


M = 


rA^ 


^eX24c2 


( 28 ) 


(29) 


we obtain 



(30) 


□ 


|W(«,1) - «l)ll < 8ce| exp] - (j^ 

This proves Theorem 1. 

We have repeatedly used the following lemma in the proof of Theorem 1. 

Lemma 2. If G*^*^(so) = 0 for some sq € [0,1] and for all k = 1... M, then 
(i) cj)^'^\so) is parallel to 4 >{sq) for all k = 0,..., M, 

(a) (so) = 0 for all k = 1,..., M, 

(in) 'Pk\so) is parallel to (j){so) for all k = 0,..., M and I = 0,..., M — k. 

Proof. Gcj) = 0, so = 0 for all k. Applying this for fc < M and evaluating the result at 

sO) the derivatives of G vanish, thus = 0 and therefore is parallel to ^ at sq, which proves (i). 

To prove (ii), use the Cauchy formula 




(31) 


where T = {z £ C | |z| = A/2} is a fixed curve around 0. Taking the /cth derivative of Eq. (31) and evaluating it at 
So, we see that the derivatives of G~^ also disappear. 

To prove (iii), we proceed by induction on k. By (i), the claim is true for k = 0. For fc > 0, we have = 

+ i{<G ~^By (i), the first term is parallel to at sq- The second term consists of derivatives of G~^ 

and derivatives of ■ By (ii), the derivatives of G~^ vanish at sq, so that the only remaining term is iG~^(f)yyy. 

But this term also vanishes at Sq by the induction hypothesis. This proves (iii). □ 

In the remainder of this section, we derive the bound on the norm of (j)M-i which was used in the proof of Theorem 1, 
following the analysis in [45]. First, we recall two technical lemmas from [45], which will be used repeatedly. 

Lemma 3. Let p, q be non-negative integers and r = p q. Then, 

^ / k'\ \ (] 4- n\U Ic — biEp-t-'rlld’ta 

(32) 


/fc\ [(Z +p)!(fc — Z + (7)!]^+“ ^ [(fc + r)!]^+“ 47r2 

^ Wy (/ + p + l)2(fc — Z + g + 1)2 ~ (fc + r + 1)2 3 


1=0 


Proof. Notice that if r = p + g, then 

^k 


[(Z+p)!(A) + g-Z)!]i+“ = 


_ (i)p+c!r: ^ p 


(k+r\ l+“ 
\l+p) 


1-)!] 


1 + Q; 


(33) 


To upper-bound the summation, divide the sum into two parts at [(fc — p + g)/2j. If I < [{k — p + g)/2j, then 
(A: —Z-l-g-l-1) > [(A:-|-p-|-g)/2j -I- 1. Otherwise, if Z > [(Zc—p-|-g)/2j, then (Z-|-p-|-l)> [(Zc-|-p-|-g)/2j-|-1. Therefore, 


E 

1=0 


^ 2 E 


(Z -l-p -I- l)2(fc - Z -I- g -I- 1)2 + 1)^ {[ik + P + g)/2j -I- 1)2 ■ 


(34) 


This can be upper-bounded by 47r2/3 as fc-|-p-|-g-|-l < 2([(fc-|-p-|-g)/2j -1-1). This finishes the proof of Lemma 3. □ 














We now use Lemma 3 to prove that if A and B are Gevrey-class, then their product is also Gevrey-class. 
Lemma 4. Let A(s), B(s) (s G [0,1]^ he smooth paths of either vectors in H or operators in B{‘H) satisfying 


9 


< Cd'^ 


[{k +p)!]^+“ 
{k+p + 1)2 


and 


||5(G|| < Ef 


[{k + g)!]^+“ 
{k + q + 1)2 


for some non-negative constants C,d,E,f, non-negative integers p^q, and for all k>0. Then, 

^2 \(i, I ^')jll+a 

r 




for all fc > 0, where g = maxid, f) and r = p-\- q. 

Proof. We have 

||(Ai3)«|!<^Q)p«||||i?('=-')||, 

so inserting the bounds (35) and upper-bounding d and / by g, we obtain 


I 


[(Z -hp)!]^+“ [{k-l-\- (?)!]^+“ 
\l) (Z+p+l)2 (k-l + q+1)^ ■ 


(35) 


(36) 


(37) 


(38) 


Using Lemma 3 to upper-bound the r.h.s. of this expression proves Lemma 4. □ 

Next we give a bound on the derivatives of the pseudo-inverse G~^. As G is non-invertible, the proof consists of 
two steps: first reducing the problem to the invertible case, then proving that the inverse of an invertible Gevrey-class 
operator is again Gevrey-class (assuming that the inverse is uniformly bounded). 

Lemma 5. If G satisfies Eq. (14), then for all k>0, 




[fc!] 


1+a 


(fc + l)2' 


Proof. First, write the pseudo-inverse using the Cauchy formula, 




(39) 


(40) 


where F = {z G <C||z| = A/2} is a fixed, s-independent curve. Taking the fcth derivative of Eq. (40) (with respect to 
s), we get 


Thus, the norm of the pseudo-inverse can be bounded by the triangle inequality, 

||(G-i)('=)||<max||[(G-z)-i]«||. (42) 

z^r 


Note that G — z is invertible and ||(G — z) ^|| < 2/A for z G F. We now show that (G — z) ^ is also Gevrey-class, 
more precisely that 


[(G-z) 


-11 (fe) 


2 

2 4tt^\ 

< — 

ztKc - 

“ A 

/A 3 j 


[fc!]i+“ 

(fc+l)2 


(43) 


for fc > 0. To show this, we proceed by induction. For fc = 0, the bound obviously holds. Taking the fcth derivative 
of (G — z)(G — z)~^ = 1, we get 


[{G-z)-^f^ 


(G - z)-i Q) (G - [(G - z)-i] 


(44) 
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Using the induction hypothesis and collecting terms (notice that I > 1 and k — I < k — 1), we get 


[{G-z) 


-11U ) 




A VA 

Using Lemma 3 to upper-bound the sum in (45), we get 

\-ii ('=) 


k\ [l\{k -1)1]^+°^ 


ij ii + i)^{k-i + iy 


(45) 


2 

f 2 

47r^ \ 

< — 

—Kc - 

- A 

V A 

3 7 


[fc!] 


1+a 


[(G-z)-i]' 

This proves (43). Substituting this bound into Eq. (42) proves Lemma 5. 

Next, we prove that the ground state is also Gevrey-class (with the special choice of the phase as above). 
Lemma 6. If G satisfies Eq. (14), then the ground state (j) satisfies 

k 


(46) 

□ 






(47) 


for all fc > 0, where K,c and a are defined in Eq. (14) and A is the minimal gap of G. 

Proof. We proceed by induction on fc. For fc = 0, (47) just reads ||(/)|| < 1. For fc > 0, notice that G(j) = 0 and 
= —G~^G^^^(f> since the phase of 4> is chosen such that (^4> (f^ = 0. Therefore, 

(fc-i) 






Expanding the derivatives, we get 


fe-i 




i=0 


fc- 1 


G-1g(i) 


H {k-l-1) 




(48) 


(49) 


The right hand side can be bounded using the induction hypothesis as the higest derivative of (j) appearing there is 
the (fc — l)th. For that, we first derive a bound on the norm of the derivatives of G“^G^^i. This can be done by 
applying Lemma 4 to G^^i and G~^ and using Lemma 5 to obtain 


jG-iG(i)](fc) 

for fc > 0. Substituting this bound into (49), we obtain 


2 471^^+^ [(fc + l)!]^+“ 
- * ^ A 3 j (fc + 2)2 


(50) 


k-1 


<E 


1=0 


k-l 


k-l\ 2 47r2y ‘ [(/c _;)!]!+“ 2 f4:Tr'^\ \ [d] 

Kc- " 


Kc- 


A 3 y {k-l + l)"^ \ A V 3 y / (Z-f 1)2' 


lH-a 


Notice that Z < fc — 1, so 




3 y V ^ y {k-i+ir + 


Thus, using Lemma 3, we obtain 




< Kc— 


[fc!]’ 


A V 3 y y (fc-f 1)2’ 


which proves Lemma 6. 


(51) 


(52) 


(53) 

□ 
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We are now in the position to bound ||^m-i||- Instead of bounding it directly, we prove a general bound on all 
The desired bound is obtained then by setting j = M — 1 and fc = 1. 


Lemma 7. For all j, k > 0, the vectors (pj and scalars ipj defined in Eq. (18) satisfy 

Jl+a 

{k + j + 1)2 

where the constants Ai, A 2 and A^ can he expressed with the constants appearing in Eq. (If): 




1+Q; 

J I - 


Ai — 




Aq — 2 c 


K ( 47r2 
A 


Ao = 4c- 


/47r 

A3 


.2\ 5 


(54) 


(55) 


Proof. We proceed by induction on j, using the recursion in relation (18). We first bound \Lpj\ using the induction 
hypothesis, then bound for fc > 0 before bounding 

Base case. We have (po{s) = 1 and <(>o(s) = (({s), so (54) holds for j = 0 since 


^ K f 47r2 

Ai > 1 and A 3 > 2c— - 

^ - A V 3 

Bound on \<Pj\, j > 1. Iv^jl can be bounded by the maximum value of the integrand in Eq. (18), 


\ f ,\ < IIG- 




(56) 


(57) 


Using the bound on ||d)|| from Lemma 6, IjG ^|| from Lemma 7 and the induction hypothesis on ||0j_i||, we get 


, |.2 1 ,-i [(j)!]^+“ ,,[(j)!]^+“ 

1^,1 - 413 - 


(58) 


since 


, .^321 

Bound on ||‘/Jj^^||. We now bound for fc > 0. First, from the induction hypothesis, 

||i(fc) II ^ A A j-l Ak+l li^ 4- j)'.] 

\\<P,. 3 \\ < a,a, a, 

Then, using Lemma 4 and Lemma 5, we get that for all fc > 0, 


(g-'0,-i) 


(fc) 


< 




nl+a 


3 A 


(fc+j + 1)2 ■ 


(59) 


(60) 


(61) 


Moreover, from Lemma 6, 


Since (pj = 


Lemma 4, Eq. (61) and (62) imply that 


1+a 


I.r ■>! = I < 


l+a 

2 )~ 


for fc > 0. Changing fc + 1 to fc gives 


L.(fc)|/ . AJ-I Aklik + + ^ AJ Aklik+jV-] 


l+a 


since 1 > 


(62) 


(63) 


(64) 
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Bound on We now bound Using the bound on from Lemma 6 and Lemma 4 implies 




1+q; 


(fc+j + 1)2 


(65) 


Finally, using Eq. (18), 


||(/)^-^^|| < 2 max 




{g 


(k) 


( 66 ) 


Hence, since 


we obtain 



and 


47r2 2 ^ 


ii-^rii 


< A.AiAl 


[(fc + j)!]i+" 

{k + j + 1)2 ’ 


(67) 


( 68 ) 


which finishes the proof of Lemma 7 and hence the proof of Theorem 1. 


□ 


Locality of local adiabatic change 

We show in this section that G„ and G„, as defined in Eq. (7) and in Eq. (12) in the main text, generate basically 
the same dynamics. The proof relies on G(0) being frustration-free, and a runtime of r = 0(log^''’“(A^/e)), because it 
turns out that the achieved locality scales linearly with the runtime. We also use the Lieb-Robinson bound [32, 49-51], 
which states that if 77 is a local (possibly time-dependent) Hamiltonian with uniformly bounded interaction strengths, 
U(t) is the unitary evolution generated by H, and Oa, Ob are operators supported on regions A, B, respectively, then 

\\[U{t)OAU'^{t),OB]\\ < cmin(|H|,|H|)||OA|j||OB|!exp( 7 t-^L), (69) 

where L is the distance between A and B, and c, 7 , v are constants depending only on the geometry of the lattice and 
the maximum interaction strength. 

The following theorem justifies the replacement of (12) with (7) in the main text, without significantly altering the 
evolution and thus the final state. 

Theorem 8. Let G(s) = “ frustration-free Hamiltonian path with 0{N) local terms such that 

I supp ^G| = 0(1), and let G be a localised version of G, i.e., 

Gis) = '^ Gf,{s), T'= |/r e T I supp ^G, supp G^) < Xrj (70) 

for some constant x o,nd adiabatic runtime r. Let ip and ip be the evolved states under G and G respectively, i.e., 

= O (^) V'(i), = O ^(f), te[ 0 ,T], •!/'(0) = V'(O) = «!>(0), (71) 

where (p{0) is the ground state o/G(0). Then, for sufficiently large x = 0(1), 

||'0(t) — iP{t)\\ < cN'^T^eO~'’^l'A'^, (72) 

where 0 , 7 , 1 ; are the constants from (69). In particular, if t = 0(log^''’“(iV/e)), then 

\\^{T)-'>P{T)\\<e/N (73) 


for sufficiently large x = 0 ( 1 ). 
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Proof. For any H C T, let Go = J2fien ^ij- Ga(t, s) be the unitary dynamics generated by Gq. Then Uq satisfies 


dtUn{t,s) = -iGn{t/T)Un{t,s), 
dsUn{t,s) = -iUQ{t,s)Gn{s/T), 

Unit, s) =T exp I y dt'Gn{t' /t) 


(74) 

(75) 

(76) 


Notice that Gt' = G and Gx = G. We write U = Uy and U = Ur- Let B be the boundary of T', that is, 
B — {^ G T \ d{X,fj.) = [x'''!} and B = T\B. Then, since G(0) is frustration-free and all terms outside of T' are 
constant, U§{t,0)cj){0) = U{t,0)cj){0) asUg = U ® and = (fiO). In other words, Gg generates the 

same dynamics as G. Thus, 

Mr) - V'(r)|| = ||G<^(0) - umw = Wm - Ugmil = ll<^(0) - U^Ugmi (77) 

where U and Ug are evaluated at (t, 0). Let V(t) = {t,0)Ug{t,0). Then, since Gg = G — Gg, Eq. (74) implies 

= iUHt,0)Gg{t/T)U{t,0)Vit). (78) 

at 

We now approximate V with a local unitary to obtain a bound for (77). Let X = {ii G T \ d{g,,B) < r} for some r 
to be specified below, and let 

Vx(t) = rexp|zy dt'uj,{t',0)Ggit'/T)Uxit',0)^ (79) 

For r = and sufficiently large x = 0(1), X and suppG are disjoint since | suppG| = 0(1), so Gxit/r) = Gx(0) for 
all t G [0,t]. Because of frustration-freeness, Gx{t/T)(j>{0) = 0, and thus the dynamics generated by Gx acts trivially 
on the initial state, i.e., Ux{t,0)4>{0) = <^(0). Thus, Vx also acts trivially on the initial state, Vx{t)(f>{0) = 4>{0). 
Hence, substituting this into Eq. (77), we get 

ii</)(o) - u^umw = \\vxm - u^ugmw = ii</>(o) - v^vmi m 


From the definition of V and of Vx , 

y (4 e) = iVliU^GgU - Ul,GgUx)V, (81) 

where Gb is evaluated at t/r. Thus, by integrating (81) and using the triangle inequality and unitary invariance of 
the operator norm, 

||</)(0)-4e</.(0)|| < / AtWU^Ggil-U]^GgUx\\= j At\\Gg - UUl^GgUxU^l (82) 

^0 Jo 

where the unitary evolutions are taken from 0 to t and is evaluated at t/r. Observe that 

a, (c/(t, s)G^(t, s)) = s)Gx{s/T)ul,it, s), (83) 

where X = T\X. Integrating (83) over s, we get 


Gg{t/T) - l}U\GB{t/T)UxU^ = -I / dstj{t,s) \g x{s/ t) s)Gg{t/T)Ux{t, s)] U^{t,s). 


Therefore, using the triangle inequality and the unitary invariance of the norm, we get 


||^/i(t) - V)(r)|| < / dt / ds U\{t,s)Gg{t/T)Ux{t,s),Gx{s/T) 
Jo Jo 

< f dt f ds 


Jo Jo 

where the second line follows from the Lieb-Robinson bound as B and X are separated by a distance r = 
proves Theorem 8. 


(84) 


(85) 

( 86 ) 

This 

□ 
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Relaxations on the assumption of a uniform gap along the path 


In this section, we show that the assumption of a spectral gap along the entire path of can be relaxed. 

Theorem 9. Suppose that G„(0) has a spectral gap of at least S > 0. Then Gn{s) has a spectral gap of at least q^S 
for all s G [0,1], where go satisfies that 1 > Q\ > gol (with Q\ as in Eg. (3) in the main text). In particular, a 
uniform gap as defined in the main text implies a constantly lower bounded gap along the entire Hamiltonian path in 
the given algorithm. 

Proof. Since G„(0) is positive semidefinite and has a non-trivial kernel, the spectral gap condition of G„(0) is equivalent 
to G„(0)2 > 5G„(0). Note that 1 > > qol (with A„^rn as in Eq. (9)). Let X„(s) = 

Then, 


Gnis )=^ n n 








( 

n 

] 

^n,L(s) 

Pu 

( 

n 

\ 

^n,L(s) 

y m^n 

) 


\ m^n 

/ 




— '?O^n.n(a)Gn(0)A„_^(s) 

= X„(s), 


(87) 


( 88 ) 

(89) 

(90) 


where we used in the second line ||>l„^ii(s)|| G (Zo ^^otice that G„(s) and X„(s) have the same kernel and are both 
positive semidefinite. Thus, the gap of G„(s) is lower bounded by the gap of X„(s). But since G„(0)^ > (5G„(0), we 
also have Xn{sf > qQSXn{s) since A~^^(s) > 1. Thus, G„(s) has a spectral gap of at least A„ > g^S. □ 


Gibbs state preparation in the non-commuting case for high temperatures 

The algorithm we presented to prepare a purifiaction of the Gibbs state of a Hamiltonian used explicitly that the 
Hamiltonian is a sum of commuting terms. Thus, one may wonder if one can apply it directly to Gibbs states of 
non-commuting Hamiltonians H. The genaral answer is no. The reason is that even though a parent Hamiltonian 
can still be defined as 


G"'(/3) = 5] G"'(/3) = ^ (91) 

MCT AtGT 

now the terms are not local (hence the superscript nl), and the norm of each term may be exponentially large in N. 
Thus, adiabatic state preparation using (91) directly takes exponential time. However, in this section we show that 
for sufficiently high, but constant temperatures, one can approximate G"* by an r = 0(logiV)-local Hamiltonian G’’ 
which is a sum of 0(poly(7V)) terms. We also show that in this case, G"* (and thus also G*") has a H(l) spectral gap 
and 0{N) norm. Because of the existence of the gap, the ground state of O'" is a good approximation of the ground 
state of G”^ 

Using the adiabatic theorem, the following algorithm runs in 0(poly(iV)) time for high enough (but H(l)) temper¬ 
atures and gives a good approximation of the purification of the Gibbs state of a non-commuting Hamiltonian: 

1. Prepare the ground state of G’'(0) = 

2. Calculate G’'(/3) 

3. Prepare adiabatically the ground state of G’’(/3) 

We first use the cluster expansion [52] to construct the approximating Hamiltonian G'’. We also show that the 
norm of G”* is 0{N). Finally, we show that the gap of G”* is U(l). For simplicity, assume that H = X)agA ^ 
of nearest-neighbour interactions, although the results and proofs generalise to other types of interactions. We also 
assume that |j/iA|| < 1. 
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Cluster expansion. We now show that C?"* can be approximated by an r = 0(logiV)-local Hamiltonian G’’. More 
precisely, we show the following result. 

Theorem 10. For sufficiently small (but constant) ffi there exists an r = 0{logN)-local Hamiltonian with 
0(poly(A^)) terms such that 

||G"'-Gni <0(l/poly(lV)). (92) 

Moreover, the terms of can he calculated in 0{poly{N)) time. 

For any function / defined on the subsets of A, define the Mobius transformations 

fiil) := ^ /(0), (93) 

ecn 

/(H) := ^(-l)l^\®l/(0). (94) 

ecn 

It is straightforward to check that the following Lemma holds [53]. 

Lemma 11 (Mobius inversion). 


/ = / = /• (95) 

For any H C A, let Hq = /m(^) = e^^^P^e^^^ for any p G T. Using Lemma 11, one 

can express G)) as 

g;:' = ffiiA) = ^ ffiiil). (96) 

ncA 

This so-called cluster expansion has many interesting properties. 

Lemma 12. Let p G T. If LI C A is sueh that p is disjoint from LI and Q C A is disjoint from LI, then 

/^(HU0)=O. (97) 

Proof. We have 

/^(HU0)= (-l)l^\^'l(-l)l®\®'le^^e'P^e-2'^^e'P^e'^^e' = 0, (98) 

n'cn, 0 'ce 

since Hqi commutes with P^ and with Hq', and the sum over H' is 0 . □ 

Lemma 12 states that is non-zero only for connected sets of edges that, in addition, contain p. Another interesting 
property of is that its norm can be bounded as follows. 

Lemma 13. For any LI G A and any edge p, 

Il/^(l^)||<(e4^-1)1^1. (99) 


Proof. Expanding the exponentials, one gets 

©CQ wi,W2,W3^Q* 


(-/3)l’"il • (2/3)l“"l • (-/3)l“3l 
jwij! • |w2|! • \w 3 \l 


h P h P h 


( 100 ) 


where 0* is the set of all finite sequences of elements of 0, and for any w G 0*, [tcj denotes the length of w and 
hw ^Ai • ■ • if rc (Ai,..., A|.j^|). 

Consider the set A = supp(wi) Usupp(u> 2 ) Usupp(u; 3 ). If A 7 ^ H, then the alternating sum in (100) over all 0 such 
that A C 0 C H is 0. Thus, 




E 

Wi,W2, W3 G LI* 

supp(wi) U supp(u> 2 ) U supp(ui 3 ) = LI 


fj\wi\ . ( 2 / 3 )I“ 2 | . p\w3\ 
jwij! • |w2|! • |w3|! 


( 101 ) 
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But this is exactly a Mobius transform, so 


\\Um< ^ 

©cr 2 Wi,W2,'W2^&* 


p\^l\ . (2/3)l’"2| . ^I«"3| 
|wi|! • |W2|! • \W3\\ 


(_l)|a\e|g4/3|e| 

ecn 


(e4/3_i)|n|. 


( 102 ) 

□ 


We are now in a position to prove Theorem 10. 

Proof of Theorem 10. Using (96), we can write Gff as a sum of local terms where the norm of the terms decay 
exponentially with their support. As the number of terms with a given size is bounded by the lattice growth constant 
[54], ||G'(]^|| =0(1) above some temperature. Indeed, let r] be the lattice growth constant, so that the number of sets 
of connected edges containing ^ and of size M is bounded by Then, 


\\G-^\\ < E ll^(^)ll < E - I)""' = OW 

OCA Af >0 


(103) 


if f3 is sufficiently small (but constant). In this case, G'ff can be approximated by an r-local operator by omitting 
all connected sets U of size at least r. The error of this approximation is 


\\g;^ - g;ii < E \\um\ < E - 1 )"" 

|f2|>r M>r 


l-ev (e4/3 - 1 ) 


1 -y’ 


(104) 


where y = — l). Therefore, above some constant temperature, the cluster expansion can be truncated at 

|fl| < O(logiV), giving an error of 0(1/ poly(lV)). This results in a 0(log(iV))-local Hamiltonian 


G"’ = E E 

M |n|< 0 (logAr) 


(105) 


with 0(poly(iV)) terms. Note that this Hamiltonian can now be calculated in 0(poly(A)) time. Indeed, there are 
0(poly(A^)) terms /^(H), and each term can be evaluated in 0(poly(iV)) time since there are at most 0(poly(iV)) 
subsets of each H. □ 

Gap of It remains to be shown that at sufficiently high (but constant) temperatures, the parent Hamiltonian 

is gapped. 

Theorem 14. For sufficiently small (but constant) (3, has a spectral gap o/H(l). 

Proof. To show the existence of a gap, we use that G"* > 0 is frustration-free, so it is enough to show that 

(Gn/)2 > ^Qnl 

for some A = H(l). Expanding G"\ (106) is equivalent to 

E(^m)' + EGm'gd'>E^Gm- (107) 

Using Eq. (104) with r = 1, we get that G)) is close to = G((=^ for high temperatures and thus it is gapped and 

the gap is close to 1. Therefore, it is enough to show that for some other constant A' < 1, 

Eg"'g"'>-Ea'g;]/ (108) 

We upper bound the r.h.s. by lower bounding G”* as 

E g”' + g/' 

U r d{n,i/)=r 


( 109 ) 
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where x = 2 e ^Cdx'^ = 0(1) is the number of times a single term is counted, and d is the dimension of the lattice. 
Therefore, it is enough to prove that for a given r and any pair fj,, v with d(/i, v) = r, 


G nl /^n 
IJ, 


> +G”'). 


( 110 ) 


Note that the kernel of the LHS of (110) is contained in the kernel of the RHS. Next, G”* +G”* can be lower bounded 
by 


GO' + G 


nl 


> - 
- 2 


5( 


1 -P, 


Ker(G]l‘+G^‘) j i 


( 111 ) 


since G))' + G))' ~ Pfj. + Pi^, which has gap 1, and at sufficiently high (but constant) temperature the difference is 
sufficiently small. 


^1021 


+ Xn and GO = 




+ x,^. 


G^j-m 


and 




To lower bound the l.h.s of (110), write G))' = 

commute as they are supported on two disjoint regions, and they are positive, thus their product is also positive. The 
norm of X^, X^, is bounded by (with y = — l) as defined in the proof of Theorem 10) 


= \\G, - < l|G^ - GlG^'^w + ||Gr^/2T - iGr-’/^Tin < 3^ 


\r/2\ 


1 - 2 /’ 

since ||G0' — — y'^ by (104), and thus g\J'G'^ > —y\'r/'^^|[\ — y'^^ so —\g\J'G'^ 

2 yPGA /(i _ yy Using that above some constant temperature ||G))|| < 2, we get that 

, 1^21 1 

!- < A' —e”’' 

1 — y 2x 


( 112 ) 


< 


G"'G0' - \gIG^']\\gIG^'^ 


< 18i 


for sufficiently small (but constant) j5. Therefore for any v pair, the following is true: 


1 1 

G;)'GD' > -18y3^[1 - - PKer(G-y+G-^)] > -A'^e-^(G;)' + Gf), 


as the kernel of G"*G"* contains the kernel of G"* + G"^ This proves Eq. (110) and thus Theorem 14. 


(113) 

(114) 
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